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It is proposed that some diffraction problems can be conveniently solved by a direct 
numerical integration of the Fresnel-Kirchhoff formula. The required properties of the dif- 
fraction screen are represented by a series of numbers which can be either regular and 
periodic, or partially random. The necessary limits and integration intervals are considered, 
and the method is found to be convenient for Fresnel diffract ion and for irregularities not 
too large compared to the wavelength. Both deep and shallow modulation can be treated. 
The accuracy of the computations is verified in a simple case of sinusoidal modulation, and 
some new results are derived for random phase screens. 



1. Introduction 

Diffraction is an important phenomenon of radio 
propagation, which is likely to be encountered in the 
many instances of irregular propagation media. In 
theoretical treatments it is usually assumed that the 
wave emerges from the irregular medium with ran- 
dom fluctuations of amplitude and/or phase along 
the wave front, and that the diffraction problem 
concerns just the evolution of this initial modulation 
as the wave propagates, no additional modulation 
being introduced further on. This leads to the con- 
cept of a thin diffraction "screen" where all fluctua- 
tions of signal are supposed to be introduced. Given 
the properties of the diffraction screen it is then 
required to find the properties of the signal over an 
observation plane some distance beyond it. The 
computed diffraction pattern may subsequently be 
compared with experimental data, in order to make 
deductions about the character of the irregular 
medium. 

Even in this simplified form the problem is not 
easy. The greatest step towards its solution was the 
development of the "angular spectrum" theory 
[Booker, Ratcliffe, and Shinn, 1950; Ratcliffe, 1956] 
which, drawing analogy with the frequency spectrum 
of a time-modulated signal, provided a valuable 
physical concept of the process and leads to direct 
answers in the more simple instances — for example, 
when the diffraction screen introduces only "shallow" 
modulation and the observation plane is far away 
from it. More difficult are the situations where the 
observation plane is relatively close to the screen— 
i.e., Fresnel diffraction — or when the screen intro- 
duces "deep" modulation, which is usually taken to 
imply phase fluctuations exceeding one radian or 
amplitude fluctuations exceeding one neper. Hewish 
[1951], Bowhill [1957, 1961], and Wagner [1962] have 
attacked aspects of these problems by the angular 
spectrum method. An alternative approach, involv- 
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ing direct application of Huygen's principle through 
the Fresnel-Kirchhoff formula, has been taken bv 
Mercier [1959, 1962]. 

By dint of mathematical sophistication both 
methods have yielded valuable results; yet simplify- 
ing assumptions, such as of a Gaussian distribution 
of fluctuations, or a small-angle approximation, 
have normally had to be invoked. Moreover, it 
does not seem to have been possible to obtain any 
numerical results for random screens whose phase 
modulation exceeds an rms value of -yj2 radians. 
The approach to be offered here, on the other hand, 
is almost entirely numerical, being based on a simple 
numerical integration of the Fresnel-Kirchhoff for- 
mula. The properties of random screens are simu- 
lated by the generation of correlated random 
numbers. While the mathematics are now ele- 
mentary, a high-speed computer and appropriate 
programming skill are required. The merits of an 
analytical solution are lost, but it seems that the 
following advantages may be gained by the numerical 
approach: 

1 . There is no need to restrict the screen properties 
for mathematical convenience. Any form of screen 
may be assumed, and the appropriate numerical 
solution will be obtained in a form suitable for direct 
comparison with experimental data. 

2. The Fresnel, or small-angle, approximations 
are avoided. 

3. Deeply modulated screens do not present any 
particular difficulty. 

Thus, the lost convenience of an analytical answer 
may conceivably be compensated by the ability 
to solve a greater range of problems. 

The present task is to explore the possibilities of 
the numerical method. We shall try to establish 
confidence in its accuracy by working out some simple 
cases whose answers are known from previous work. 
We will examine the limits within which the method 
looks practicable (these arc in effect limits of com- 
puting speed) and compute some more difficult in- 
stances as examples. 
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2. Formulation 

The starting point is the Fresnel-Kirchhoff diffrac- 
tion formula 3 [Born and Wolf, 1959], which says that 
the field at B due to a source S is 



U 

where 



2\Ja 



exp jk(r-\-s) 



(cos x r — cos x s )dA (1) 



and 



dA is an element of a surface A enclosing R or 

S (or an infinite plane between them), 
r and s are the distances of R and # from dA, 
X r and x s are the angles between the normal 
to dA and the directions to R and S. 



2tt 



k=— > for wavelength X. 

A 



If the source is at infinity, 



u= -iL 



exp $r 



(1+cos x r )dA. 



(2) 



In the present work, A will be a plane distant h 
from R, as in figure 1, and the simplifying assumption 
will be made that irregularities exist in the x-direc- 
tion only, being infinitely elongated in the ^-direction. 
Diffraction in the ^-direction does not then concern 
us, and the formulation can be reduced to two 
dimensions by applying to the z -direction the usual 
Fresnel approximations: 



SOURCE 




Figure 1. Geometry of the model. 



3 Although not rigorous, the Fresnel-Kirchhoff formula is quite the most tracta- 
ble expression of Huygen's principle. The usual objection that it neglects 
boundary effects at an opaque screen will not apply here because no discontinui- 
ties will be involved, and curvature of the constant-phase front will always be 
greater than a wavelength. The derivation of (1) from the more basic Helmholtz- 
Kirchhoff equation assumes r, s, »X/2tt. 



r ~*( 1+ 5?) 



in the exponent, 



and 



q in the denominator, 
x r ~x. 
Whence, for the model in figure 2, 

exp jkq (1 + cos x) 



U- 



vix: 



2* 



dx. 



(3) 



In contrast to previous analytical solutions, no 
approximation has been made in the x-direction 
where the irregularities are. 

It is now supposed that the constant-phase front 
of a wave which originated at infinity has been 
delayed during propagation to the plane A by 
amounts Ah(x), as in figure 2. Thus, assuming 



U 



Vil 



+ 00 exp jkUh 2 +x 2 +Ah cos x) (1+cos x) 
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dx. 

(4) 



Keplacing this by a summation, and expressing 
distances in terms of the wavelength, gives 

= S SS ex P 2xi [HSe +AH/Se] (5) 



\j]H^ L 2Se 3 



where 



Al=^, AH- 

A 



=— ,H=-, and Se-- 

A A 



-Sec 
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The limits are now ±X L instead of + <» . For the 
other case of practical interest, where both S and 
R are distant h from the screen, similar reasoning 
gives 

Ai 



AX +x L exp27rj^2HSe+~j 



Se 5 ' 2 



(6) 



In general it will be convenient to express the field 
as an amplitude and a phase term: 



U=M exp (j(p). 



(7) 



Our general procedure will be as follows. The 
modulation function Ah(x) is generated, and the 
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OBSERVING PLANE, A 



Figuee 2. Phase irregularities. 
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field U is determined from (5) or (6). The modula- 
tion function is then translated along the ^-direction 
in steps of Ax' , the corresponding; U(x') being 
evaluated at each step. Thus the diffracted field 
along the observing plane will be derived, given an 
assumed initial modulation at plane A. In the 
present work only sinusoidal or random screens 
will be used, and the diffracted fields will be derived 
from (5), corresponding to the source at infinity. 

3. Preliminary Tests: Unperturbed Wave 
Front 

Before (5) can be used in computing diffraction 
problems, it is necessary to find the smallest summa- 
tion limit, X L , and the largest increment, AX, that 
can be used without undue error. They in turn 

(2X 4-1) 
specify how many terms, - — -^ — > must be taken, 

and so indicate the feasibility of obtaining solutions 
with a given computer. For this reason initial com- 
putations were made assuming the wavefront to be 
unperturbed — i.e., All=0 — in which case the result 
should be M=l, <p=0. 

The progressive summation of (5) may be regarded 
as the construction of a spiral. When X is large 
relative to H the terms become small and equal to 

AX 

in magnitude. Provided AA" is small, one 



2^X 

revolution of the spiral is completed every time X 

increases by unity, i.e., every -pr? terms, and the 



AX 



spiral has radius 



1 



Thus, if (5) is summed to 

4:T-yJX 

X L is large relative to H, 



±X L , where 

residual error with magnitude 



there is a 



1 



2t^X l 



with respect to the sum to infinity. In the more 
general case when X L is not necessarily large com- 
pared to H and the individual terms are of magnitude 

AX (1+Se) 



2V# Se 3/2 
the residual error may be obtained from figure 3 
— as a function of -77 • This was 
computed for H=5 and shows, for instance, 



that e = 2.6 percent if 



H 



9 and H=5. In this 
91 



instance it would be necessary to include "Z~ terms. 

J AX 

Fewer terms would be needed in the summation of 

(6) to the same accuracy because the individual 

terms, again at large X, fall off in magnitude as 

1/X 5/2 rather than 1/X 1 ' 2 . 

The procedure adopted for computation with (5) 

is to evaluate the series up to some X L for which e is 




Figure 3. eVH/5 in terms of X L /H. 

suitably small (usually 2 to 3%) and then to add on a 
computed error term of appropriate amplitude and 
phase. Test computations with II between 2\ and 
20, taking e=2.5 percent, gave remaining errors of 
less than 0.2 percent after the addition of the error 
term. 

In the course of the foregoing tests it was found 
that the increment, AX, had no significant effect on 
the results provided AX<}{. It seems that a rela- 
tively large interval will suffice for an unmodulated 
wavefront. However, this may not hold for modu- 
lated waves, where it will obviously be necessary to 
delineate the perturbations in sufficient detail, and 
the question of increment will be considered again 
later. 

4. Sinusoidal Screens 
4.1. Shallow Modulation 



Diffraction by a sinusoidal phase screen can readily 
be evaluated by the angular spectrum method 
[Hewish, 1951; Ratcliffe, 1956] when the modulation 
is shallow — i.e., when the fluctuations are small 
relative to 1 radian — and only first-order side- 
waves need to be taken into account. It is there- 
fore a suitable test of the numerical method. The 
modulation is of the general form 



AH=AlL cos 



2ttX 
P ' 



This calculation has been carried out by both the 
angular-spectrum and the numerical methods for 
a modulation depth (Aflo) of 0.01 wavelength and 
for different values of the period (P), the distance 
from diffraction to observation plane (H) being taken 
as 5 wavelengths. As is well known, the diffraction 
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pattern is also sinusoidal, with the same period as 
the irregularities, but the actual magnitudes of the 

phase and amplitude fluctuations (<p and M) depend 

on both P and H. <p and M are shown as a function 
of P, for H=5, in figure 4. The two methods, in 
fact, give virtually the same result, differences being 
no greater than 0.001. X L /H=9 was used for the 
numerical method, and the same results were ob- 

The " total fluctua- 



tained for AX=%, He 



and % 2 - 



tion," which is V^+^yr 2 , is independent of P when 
P is large compared to one wavelength, and figure 
4 illustrates how the fluctuation alternates between 
amplitude and phase, being mainly in the phase 
for periods ^ornewhat greater than the first Fresnel 

zone (P>V#). 

The close agreement between the two methods in 
the case of a shallow sinusoidal phase screen shows 
that the numerical method is capable of good 
accuracy in a simple instance. We shall now 
evaluate some more complex examples, for which 
analytical results are less readily available. 

4.2. Deep Modulation 

As long as the fluctuations in the screen are small 
compared with 1 radian (AH <^<^l/2ir) the form of 
the diffraction pattern is independent of the initial 

fluctuation depth, and the values of M and <p are 
simply proportional to it. When the initial modula- 
tion approaches 1 radian in depth these simple 
truths fail and then the modulation can no longer 
be considered "shallow." The results of a series 
of computations for AH =Q.l are plotted, suitably 
scaled, on figure 4, and a considerable deviation 
from the form of the shallow modulation can be 
seen. In terms of the angular-spectrum concept, 
these changes with increasing depth of modulation 
are due to contributions from side-waves of higher 



0.08 

0.07 

0.06 

0.05 

0.04 

0.03 

0.02 

0.01 



-0.01 

-0.02 

-0.03 

-0.04 

-0.05 

-0.06 

-0.07 

-0.08 




_L 



I 



3 



5 



7 



PERIOD (P) 
Figuee 4. Shallow sinusoidal phase screen. 



order that introduce higher spacial harmonics into 
the diffraction pattern. Examples with AH near 
or exceeding 0.1 — a phase fluctuation of 7r/5 radian — 
plainly cannot be treated as shallow modulations. 
In the angular-spectrum method, problems invok- 
ing deep modulation require the addition of several 
side-waves whose amplitudes are specified by a 
Bessel function. The number of required terms 
increases rapidly with increasing modulation depth. 
By this technique Hewish [1951] computed diffrac- 
tion patterns from initial phase modulations of up 
to 8 radians. In the present numerical method 
deeply modulated screens can be treated as easily 
as shallow ones, apart from the need to verify that 
the summation interval (AX) is adequately small. 
Figures 5 and 6 show the diffraction patterns com- 
puted for the conditions H=5, P=3%, and H=20, 
P=2, with initial modulation depths (AH ) up to 
four wavelengths (25 radians). These show plainly 
the distortion of the sinusoidal form by increased 
modulation depth, with periodicities down to half 
a wavelength in evidence. In general the overall 
amplitude fluctuation of the diffracted signal does 
not change greatly beyond AH =0.1; whereas the 
overall phase fluctuation (which has been plotted 
as a fraction of the input modulation <p s =2irAH ) 
tends to increase at least as far as AH =1.Q. The 
curves are similar in general appearance to those 
given by Hewish [1951], but a detailed comparison 
is not possible because he does not state his value 
of P. 

4.3. Summation Interval for Deep Modulation 

The foregoing computations were made with a 
summation interval AX=}(e. To investigate the 
adequacy of this interval, some computations were 
repeated using other intervals, in the range from 
%2 to unity. Figure 7 shows the patterns obtained 
in the case H=20, P=2, AH =0.1. The results 
are identical from intervals }U to }{, and even interval 
K incurs very little error in spite of the fact that it 
barely specifies the waveform adequately. Compu- 
tations for H=5, P=3%j Afl =l, the most severe 
case of figure 5, using AX=% 2 , }U, % and %, differed 
from each other by less than 0.0024 in either ampli- 
tude or phase; relative to a total fluctuation of 1.3 
in amplitude and of 5.8 in phase this is a negligible 
error. An even more severe case in figure 6, H=20, 
P=2, Ai?o=4, was computed using several intervals 
between %2 and %. Here the results for % 2 and }U 
agree, yet those for }{% and % were utterly different, 
as figure 8 shows. It is curious that the error should 
increase so drastically for intervals which are only 
slightly too large. Figure 9 clearly illustrates the 
rapid deterioration at two particular points of the 
pattern: intervals less than }{*> are accurate, Xs 
shows a small error, and intervals larger than }{ 5 
cause unacceptably large errors. 

The sharpness of the transition between satis- 
factory and unsatisfactory intervals has one useful 
aspect. It is possible to test an interval by simply 
repeating the computation with the interval doubled. 
It the two results agree it is highly probable that 
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Figure 7. Effect of summation interval (AX) for AH o = 0.10. 
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Figure 9. Effect of summation interval (AX) on computed 
amplitude. 

both intervals are satisfactory. According to these 
tests the interval can be so large that the phase 
fluctuation in the diffraction screen is permitted to 
change by several radians between sample points. 
This, and the need to specify the shape of the 
modulation in sufficient detail, are the criteria to 
be kept in mind when deciding what summation 
interval should be used in a particular instance. 

4.4. Discussion 

The numerical method appears to be accurate in 
computing diffraction from both shallow and deep 
sinusoidal phase screens. It is estimated that on 
an IBM 7090 computer the calculation for H=20 
using a summation interval of }U would take 0.21 
min; so even deep-modulation problems can be 



evaluated quickly and cheaply. The diffraction 
screens which arise in nature may not normally be 
sinusoidal in form; nevertheless, as an approximation 
it is sometimes convenient to consider them so, and 
t he methods outlined so far may be useful in this way. 

5. Random Screens 

5.1. Simulation of Random Screens 

The diffractio nscreens appropriate to the natural 
phenomena of radio propagation are usually random 
in the sense that it is not possible to predict their 
space and time characteristics in detail. It is usually 
assumed, though, that their statistical properties 
may be exactly defined by parameters such as the 
standard deviation and the autocorrelation func- 
tion — that is, that the screens are statistically sta- 
tionary. We will specify a random phase screen by 



r(Atf) = V(A#) 2 



,(Afl> 



_(AH) n (AH) n+m 



\AH) 



<t being the standard deviation and p m the correlation 
coefficient with respect to interval m. It is now 
necessary to generate a set of numbers having these 
statistical properties, which will represent the varia- 
tion of AH across the screen. 

A series of correlated random numbers, a, may 
conveniently be generated from a series of un- 
correlated random numbers, b, by 



dn + l^fan + gb, 



(8) 



where / and g are coefficients. If the standard devia- 
tion of the a's is to equal that of the b's, 



<?=Vi-/ 2 . 



(9) 



Then the mth correlation coefficient is given by 

Pm=j m . (10) 

This simple form will be used in the present work; 
others could be devised if required. Figure 10 
shows the values of p m computed from 500 or 1000 
numbers generated by this series for three values 
of/. They are seen to agree closely with/* 1 . 

5.2. Procedure 

Since the numbers are random it is necessary to 
use enough of them when evaluating their statistical 
properties. As a guide we note (as can easily be 
proved) that if the correlation coefficient is of ex- 
ponential form, 

p m =exp (—m/m Q ), 

then v points of the correlated series contain the same 
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Figure 10. Autocorrelation of generated random screens. 

information as r/2m independent points. For smal 
values of p, its standard deviation when calculated 
from N independent points is 



<Kp) = 



1 



VA^-3 



We therefore expect the accuracy of an estimate of p 



to be of the order of 



\2m 



where r is the num- 



ber of correlated points used to evaluate p, and rn 
is the interval at which p has fallen to about 
l/<?( = 0.37). Table 1 shows m and <r(p) appropriate 
to 1000 points, for several values of/ to be used in 
computations. The value of Q, defined as the 
distance in wavelengths at which p falls to 0.5, is 
also given assuming AX=l/8. 

Table I. Properties of correlated random numbers 



f 


Wo 


o-(p) (for r = 1000) 


0(UAX=1/S) 


0.940 


16.1 


0.19 


1.44 


.884 


8.1 


.13 


0.728 


.794 


4.3 


.09 


.390 


.707 


2.9 


.08 


.258 


.500 


1.4 


.05 


.125 



In the following computations a series repre- 
senting AH is generated from (8), appropriate 
values for/ and a(b) having been chosen in advance. 
The properties a (AH) and p(AH) are computed 
for this series. The series AH is then supposed to 
represent the phase fluctuation sampled at intervals 
AX along the diffraction plane, and the amplitude 
M of the diffracted wave is computed at intervals 
AX' (where AX' can be equal to AA r , 2 AX, etc.) 
along the observation plane. The computation is 
restricted to amplitude for the sake of brevity, this 
being the quantity most often observed experi- 
mentally. M is evaluated at 500 or 1000 intervals 



and, finally, the statistical properties a(M) and 
p(M) are computed for comparison with a (AH) 
and p(AH). The distance in wavelengths at which 
p (M) falls to 0.5 is defined as q. 

It is verified that the interval AA" is sufficiently 
small by repeating the computation with every other 
point of the series AH eliminated and with AX 
doubled. With satisfactory intervals there is no 
significant difference between the two computations. 
The intervals AX= 1/8 and AX= 1/4 were always found 
to be adequate with shallow modulation, as would 
indeed be anticipated from the foregoing sinewave 
results. 

5.3. Shallow Modulation 

To illustrate the application of the numerical 
method to shallow random screens, some problems 
will be explored which lie on the borders of approxi- 
mations commonly applied in radio diffraction theory. 
It is generally accepted, for instance, that irregu- 
larities much less than a wavelength across (#<<0) 
do not propagate and therefore do not appear in the 
diffracted signal. For irregularities larger than a 
wavelength it is well known that the statistical form 
of the diffracted signal becomes identical to that of 
the diffraction screen when the first Fresnel zone 
over the diffraction screen is so large that it includes 

many irregularities, i.e., l<<$<<Vff- It is not 
always evident, however, at what point such approxi- 
mations become valid. As conditions like these have 
no special significance in the present numerical 
method, but also present no particular difficulty, we 
shall use it to investigate 

(i) diffraction from irregularities about one wave- 
length across: #^l;and 

(ii) the transition from the Fresnel to the Fraun- 
hofer region: Q^^H. 

A range of examples has been computed with Q 
between and 2.7 and with H between 3 and 30 
(i.e., ^/H between 1.7 and 5.5). Typical results are 
given in figures 11 and 12 which show how the dif- 
fracted pattern changes with the distance (H) 
between the diffraction screen and the observation 
plane for two values of Q. The diffracted amplitude 
pattern is obviously narrower near the screen, and 
the amount of fluctuation (a (M)) is smaller. Also, 
the smaller Q, the less is the dependence of q on H. 
In figure 12, where (?<0> the correlation function 
has a different shape at the two planes, which is 
consistent with the evanescence of irregularities 
much smaller than a wavelength. Nevertheless a 
considerable depth of fluctuation remains even at 
H=30. Figure 13 shows the dependence of the 
pattern width in the observation plane (q) on that 
in the diffraction plane (Q) at a constant H=9. In 
all cases ^H^Q. It is seen that q<CQ for the larger 
values but that q^>Q for the smaller ones. Overall, 
q varies much less than Q. The depth of fluctuation 

( expressed as / A r/x ) passes through a maximum 

near Q==l. The loss of fluctuation for (?>1 is prob- 
ably because Q is not sufficiently small compared 
with the first Fresnel zone, and the loss for $<1 is 
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Figure 11. Diffraction patterns for Q = 2. -•>.;. 
probably caused by evanescence of the finer compo- 
nents. Nowhere in these computations does / T i r 

a (AH) 

reach the theoretical limit V^r (=4.44) which is an- 
ticipated in the condition 1<^<CQ<C<^H. It is 
clear that situations up to ^H=4Q remain well 
within the Fresnel region, with o-(M)<^J~2to-(AH) 
and q<CQ. 

A particularly interesting case is that for Q=0. 
In contrast to previous work [Bowhill, 1957; Rat- 
cliffe, 1956], the present calculations give g>0, 
o-(M)>0. This case has been computed for a range 
of H between 3 and 48, the calculations being per- 
formed with two different sets of random numbers. 
The results are given in table 2. There is no ap- 
parent trend with H, differences being attributed to 
statistical error. The mean of these curves is plotted 
in figure 14. The existence at Q=0 of a finite q of 
about X/4, and of a residual <r(M) o about 2a(AH), 
is probably due to contributions from signals scat- 
tered obliquely; these will obviously be correlated 
for displacements X f up to about X/4. Theoretical 
treatments that presume a small-angle approxima- 
tion will miss these contributions. 
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Figure 12. Diffraction patterns for Q = 0.62S. 



Table 2. Computations for Q = 



Random numbers 


H 


Q 


cr(M)MAH) 


First set 


{ ^ 

6 

I 12 

24 

[48 

3 

6 

\\2 

24 

148 


0. 2812 
.2687 
. 2875 
. 3000 
. 3500 
.2875 
. 2687 
.2750 
. 3375 
.3062 
. 2962 
. 0265 


1.813 
1.871 
1.827 


Second set— . 


1.719 

2.139 
2. 015 
1.771 
1.940 


Mean _. _.. _._ 


1.981 
1.943 
1. 902 


Standard deviation- ___ 




0.120 









As a summary of the behavior over the full range 
of computations, figures 15, 16 and 17 show the 
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yo(AH) 0.5 




/>(M) 0.5 




F = 0.940 

0.940 

(WITH Ax =1/4) 

0.884 



12 3 4 5 6 7 

x' (UNITS OF A/8) 
Figure 14. Mean of 10 computations for Q = 0. 
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Figure 13. Diffraction patterns for 11 = 9. 

694 




0| I I I I I I I I I I I I I ! I I I I I I I II I I 



q 0.5 




■ MEAN OF Q= (VARIOUS H) 



1 I I I I 1 I I I I 1 I 1 I I 1 1 I I I 1 I I 



2.0 



3.0 



q 0.5 





1 ! 


1 1 ! ! 


1 1 


1 


1 

Q = 2.5 


1 1 
-27 








So 


















Q= 1.4 


-1.6 






^z 


- jt + — - 


~^+ 


Q = 0.625 

= 0.31-037 


— 




• Q= 0.238 








— 


— 




• Q = O.II2 






Q--0 












. 1 1 


1 1 


I I 1 1 


1 1 




1 


1 1 


1 


2 
1 


3 4 

I ! 1 1 


5 

i 1 


1 


6 


7 







3 



6 9 12 15 20 25 30 

H 

Figure 15. Shallowly modulated random screen: width of 
diffracted amplitude pattern. 

4 



3 — 



CfJM) 

cr(AH) 








o Q = 2.5-2.7 

X Q = 1.4- 1.6 

WAGNER (1962) 

■ BOWHILL (1957) 



12 3 4 

yR"/Q 

Figure 16. Change of amplitude fluctuation with yH/<f>. 
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(b) VARIATION WITH (q-0.30)/Q 



COMPUTED • 

THEORETICAL O 



0.5 1.0 

(q-0.30)/Q 
Figure 17. <r(M) /a (AH) for shallow modulation. 

width of the amplitude pattern (q) and the depth of 
amplitude fluctuation O(M)). Figure 15 illustrates 
the behavior of q: 

(i) If Q <1, q increases with increase of Q. How- 
ever, it changes relatively slightly with H, passing, 
in fact, through a maximum. 

(ii) If Q >1, q changes more slowly with Q but 
depends more strongly on //. It appears that 
q oc^H, at least to Jh=2Q, when Q=2.5. 
The fluctuation depths computed for Q >1 are given 
in figure 16 as a function of^/H/Q. a(M) is expected to 
fall off at the smaller values of ^H/Q because ir- 
regularities very large compared to the early Fresnel 
zones will make negligible contribution. By similar 
reasoning, a(M) must be zero when V^/Q is zero. 
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The slope found here is very different from the 
theoretical results of Bowhill [1957] and Wagner [1962] 
which are also plotted. It must be realized, though, 
that the three examples are for different forms of 
irregularity: BowhilFs for a Gaussian correlogram 
falling to 0.61 at Q; Wagner's for a "medium cut- 
off" power spectrum; and the numerical method for 
a correlogram of form/™ falling to 0.50 at Q. 

It was not possible to find any simple empirical 
relation that would describe the variation of a(M) 
in terms of H and Q over the whole range of Q. A 
relation between a(M)/a(AH) and g derived from the 
same computation is apparent, however. As indi- 
cated in figure 17a, the points derived for a given Q 
tend to lie on a straight line of slope proportional to 
1/Q which passes near the point computed for Q=0. 
This rather curious result is further illustrated in 
figure 17b, which shows how a(M)/a(AH) is simply 
related to a quantity (q — 0.30)/Q. The scatter is 
probably no more than statistical error, the points 
which deviate most having been derived from small 
values of Q. For diffraction in the Fraunhofer 
region (^H >>(?) it is to be expected that 
a(M)/a(AH)=4A4, and that q = Q (both being large 
compared to unity) so that (q — 0.30)/Q=l. The 
trend of figure 17b indeed passes through this point. 

In illustrating the use of the numerical method 
with shallow random screens, the following two 
results have emerged: 




M cr(AH) 



Figure 18. Deep random modulation. 



(i) Even very small irregularities produce a diffrac- 
tion pattern with finite fluctuation depth and finite 
pattern width. 

(ii) Fraunhofer diffraction — in which the forms of 
the diffraction screen and of the amplitude pattern 
in the observing plane are identical, and in which 
the relative fluctuation depths are simply related — 
cannot be considered to apply until the first Fresnel 
zone is many times the pattern width. 
Such results will have a bearing on the interpretation 
of experimental data of LF and VLF fading, since a 
100 kc/s signal received at the ground is 30 wave- 
lengths from a diffraction screen at a height of 90 
km and a 15 kc/s signal is a mere 3.5 wavelengths 
from a screen at 70 km. 

5.4. Deep Modulation 

To illustrate the method's application to deep 
random phase screens, one example has been com- 
puted which illustrates how the pattern width and 
depth of amplitude fluctuation change with increasing 
depth of phase modulation. It is taken that the 
diffraction and observation planes are nine wave- 
lengths apart, and the phase perturbation is of small 
horizontal scale: Q=0.369. The results are shown 
in figure 18. The effect of increasing a (AH) is ini- 
tially to reduce the width of the amplitude pattern, 
to reduce the mean amplitude, and to increase the 
fluctuation depth. Beyond <t(AH) = 1, however, 
these quantities are virtually independent of modula- 
tion depth. Presumably, additional contributions 
to the fading now go into components with structure 
too small for them to propagate a distance of nine 
wavelengths. An additional point of interest is that 
a(M) =M 2 when a (AH) is large. The interval AX= 
% was found to be adequately small in all these 
computations. 

5.5. Discussion 

Because of the need to obtain an adequate statisti- 
cal sample, computation with random screens is 
much lengthier than it is with sinusoidal ones. For 
instance, the computation takes 11.6 min in the 
case H=20, AX=%, when 1000 amplitude points 
are taken. A limited amount of random screen 
computation is probably justified, however, because 
of the closer approach to experimental reality. A 
set of computations, such as those for shallow modu- 
lation, can clearly be applied to practical propagation 
studies. If it seems reasonable to assume an initial 
modulation of phase only, and if the present form 
of correlation coefficient is taken, then figure 15 
relates g, Q, and H. An observation of q and an 
estimate of H then lead to Q. Having measured 
a(M) it is then possible, from figure 17, to determine 
a (AH). It may also be possible to find a (AH) from 
a(M) and a(<p)- — though this has not been done here. 
In that event one could work back to a determination 
oiH. 

A number of improvements and extensions are 
possible. Computational refinement may speed up 
the computer program (whose block diagram is given 
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in appendix 1). It would be logical to make com- 
putations for different forms of correlation coef- 
ficient, for a point source, and for oblique-incidence 
propagation. The model could conceivably be 
extended to two-dimensional (or even three-di- 
mensional) irregularities. Finally, it may be advan- 
tageous to evaluate the phase as well as the amplitude 
properties of the diffracted waves. 

6. Conclusions 

The object of this work has been to investigate a 
proposed numerical technique for solving diffraction 
problems. The technique in its present form appears 
to be appropriate for radio problems where the 
irregularities in the diffraction screen are a few 
wavelengths across and where the observation plane 
is not more than a few tens of wavelengths away. 
These are conditions of practical concern in the 
ionospheric propagation of LF and VLF signals. 
It has been shown that the technique gives accurate 
results with sinewave modulation across the screen, 
and it has been found that, even for deep modulation, 
integration intervals can be larger than might be 
thought at first sight. It has been possible to derive 
original results for random screens observed in the 
Fresnel region. 

Although the numerical method is not always the 
best one to use, and is limited by considerations of 
computing time, there are cases where it is very 
convenient and it may be particularly valuable in 
the interpretation of experimental diffraction data. 
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8. Appendix 1 . Block diagram of 
computer program 

Road basic parameters 
H, (1 > 1/AX,< 2 > XlIH,™ aX'JaX,m F (5) . 
Required standard deviation of AH^ ) = a T (AH). 
Number of values of M (7) to be computed =N M . 
Number of autocorrelation coefficients, p, to be 

evaluated for All and for M, and interval 

between successive values. 
Whether or not every other value of All is to be 

omitted. (8) 
Starting parameter For random number generator. 

Compute 
AX/H,AX^H. 

N^ = (l+X L /AX+l/2AX) 
A Y f 

N$o>=(2N-l)+^(N M -l) 



G=T/1-F 2 



Generate random screen 



If all AH used 



I 



Generate random number (B) with 
standard deviation a r (AH) 
Compute AH n+1 =FAH n +GB 



repeat for (N T +50) 
Store AH except for first 50 values (11) 



If every other AH used 

i 

Generate B x 

Compute D=FAH n +GB l 

Generate B 2 

Compute AH n+1 = FD + GB 2 



repeat for 



N T +50 



Store AH except for first 25 points 



(id 
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Find statistical properties of AH 
Compute mean, AH 

standard deviation, cr(AH) 
autocorrelation coefficients, p(AH). 
Output these. 



Compute 

1/Se {12 \ (HSe)/2S<?' 2 , HSe, for ^=0 (~\ ^±^ 

Generate 

sine and cosine tables 0(3°)357° and store 
Compute 



tfi=S ^S cos 2x (HSe+^y C7 2 =S |^8 sin 2n (HSe+^ 



for 



X_-X L fAX\ +X> 
H H 



(#> 



# 






AX 



cos 



2* (HSe+^y 5C7 2 =1] ^g sin 2* (HSe+^) 

X X l (AX\Xl+\ 
0V H~ H \H ) H 



M=^^(u l +du 1 r+(u 2 +5u 2 r 



Shift 



AX' 
values of AH by -r-^- intervals. 



Repeat for iVjw 



Find statistical properties of M 
Compute mean, M, 

standard deviation, <r(M) 
autocorrelation coefficients, p(M). 

Output these 

— Output next random number parameter. (14) 



1 H is the distance in wavelengths between diffraction and observation planes. 

2 AX is the distance in wavelengths between sample points in the diffraction 
plane. 

3 Xl is the limit of summation in the diffraction plane. 

4 A A"' is the distance in wavelengths between points in the observation plane 
at which results are to be computed. 

5 F is the parameter determining the width of the correlation pattern at the 
diffraction screen. 

fl AH is the phase deviation. 

7 Mis the amplitude at the observation plane. 

8 Alternate values of AH can be omitted to test whether AA" is sufficiently 
small. 

9 N c is the number of values of \X\ required. 

10 Nt is the total number of diffraction screen samples required. 



u Early points are dropped to ensure statistical uniformity in the series. 

12 S e= Vl+fl 

1 3 5 Ui, 8U2 are the correction terms for finite Xl. 

14 The next random number parameter is saved so that computations can be 
repeated with an independent set of numbers. 

Note: The program is in Fortran. Tape storage is used for extra capacity. 
Assuming Nm = 1000, N c can be up to 4000. 



(Paper 67D6-298) 
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